Nonlinear properties of left-handed metamaterials 
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We analyze nonlinear properties of microstructured materials with negative refraction, the so- 
called left-handed metamaterials. We consider a two-dimensional periodic structure created by arrays 
of wires and split-ring resonators embedded into a nonlinear dielectric, and calculate the effective 
nonlinear electric permittivity and magnetic permeability. We demonstrate that the hysteresis-type 
dependence of the magnetic permeability on the field intensity allows changing the material from 
left- to right-handed and back. These effects can be treated as the second-order phase transitions in 
the transmission properties induced by the variation of an external field. 
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Recent theoretical studies 0, 0, and experimental 
results ^ in 1^ have shown the possibility to create novel 
types of microstructured materials which demonstrate 
the property of negative refraction. In particular, the 
composite materials created by arrays of wires and split- 
ring resonators were shown to possess the negative real 
part of the magnetic permeability and dielectric permit- 
tivity for microwaves. These materials are often referred 
to as left-handed materials or materials with negative re- 
fraction. Properties of such materials were analyzed the- 
oretically by Veselago long time ago , but only recently 
they were demonstrated experimentally. As was shown 
by Veselago , the left-handed materials possess a num- 
ber of peculiar properties, including negative refraction 
for the interface scattering, inverse light pressure, reverse 
Doppler and Vavilov-Cherenkov effects. 

So far, all properties of left-handed materials were 
studied in the linear regime of wave propagation when 
both magnetic permeability and dielectric permittivity 
of the material are assumed to be independent on the 
intensity of the electromagnetic field. However, any fu- 
ture effort in creating tunable structures where the field 
intensity changes the transmission properties of the com- 
posite structure would require the knowledge of nonlinear 
properties of such metamaterials which may be quite un- 
usual. In this Letter we analyze, for the first time to our 
knowledge, nonlinear properties of left-handed metama- 
terials for the example of a lattice of the split-ring res- 
onators and wires with a nonlinear dielectric. We show 
that the effective magnetic permeability depends on the 
intensity of the macroscopic magnetic field in a nontrivial 
way, allowing switching between the left- and right-handed 
materials by varying the field intensity. We believe that 
our findings may stimulate the future experiments in this 
field, as well as the studies of nonlinear effects in photonic 
crystals where the phenomenon of negative refraction is 
analyzed now very intensively 8, 9j . 

We consider a two-dimensional composite structure 
consisting of a square lattice of the periodic arrays of 
conducting wires and split-ring resonators (SRR), shown 
schematically in Fig.^ We assume that the unit-cell size 
d of the structure is much smaller then the wavelength of 



the propagating electromagnetic field and, for simplicity, 
we choose the single-ring geometry of a lattice of cylin- 
drical SRRs. The results obtained for this case are qual- 
itatively similar to more involved cases of double SRRs 
properties. This type of the microstructured materials 
has recently been suggested and built in order to create 
left-handed metamaterials with negative refraction in the 
microwave region Q . 

Negative real part of the effective dielectric permit- 
tivity of such a composite structure appears due to the 
metallic wires whereas a negative sign of the magnetic 
permeability becomes possible due to the SRR lattice. 
As a result, these materials demonstrate the proper- 
ties of negative refraction in the finite frequency band, 
loq < Lo < min(a;p, Ci;||„i), where wq is the eigenfrequency 
of SRRs, cl;||„i is the frequency of the longitudinal mag- 
netic plasmon, ojp is an effective plasma frequency, and 
w is an angular frequency of the propagating electromag- 
netic wave, {£, H) ^ (E, H) exp (iujt). The split-ring res- 
onator can be described as an effective LC oscillator (see 
Ref. ■^ith capacitance of the SRR gap, as well as 

effective inductance, and resistance (see the upper insert 
in Fig.P). 

Nonlinear response of such a composite structure can 
be characterized by two different contributions. The first 
one is an intensity-dependent part of the effective dielec- 
tric permittivity of the infilling dielectric. For simplicity, 
we assume that the metallic structure is embedded into a 
nonlinear dielectric with permittivity that depends on the 
intensity of electric field in a general form, ejj — e£)(|E|-^). 
For detailed calculations presented below, we take the 
linear dependence that corresponds to the Kerr nonlin- 
earity. 

The second contribution into the nonlinear properties 
of the composite material comes from the lattice of res- 
onators, since the SRR capacitance (and, therefore, the 
SRR eigenfrequency) depends on the strength of a lo- 
cal electric field in a narrow slot. Additionally, we can 
expect a nonlinear eigenfrequency detuning due to a reso- 
nant growth of the charge density at the edges of the SRR 
gap. The intensity of the local electric field in the SRR 
gap depends on the electromotive force in the resonator 
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FIG. 1: Schematic of the composite metamaterial structure. 
Lower insert shows a unit cell of the periodic structure while 
the upper inset shows the SRR equivalent oscillator with the 
parameters used in the derivation. 



loop, which is induced by the magnetic field. Therefore, 
the effective magnetic permeability /^off should depend 
on the macroscopic (average) magnetic field H. 

For the polarization shown in Fig. ^ (lower insert), the 
main contribution to the dielectric function is given by 
the array of wires. When the wire length is large enough, 
so that the frequency of the fundamental (dipole) mode 
of an individual wire becomes much less than uj, only the 
resistance and inductance of the wires give a contribution 
into the lattice impedance. Hence, the Ohm law for the 
current can be written in the form, 



1 + iuaSL, 



-E', 



(1) 



where is the electric current density in the wire, E' is 
the local electric field, a is the conductivity of the wire 
metal, ~ 2c~'^ \n {d / r) {d ^ r) is the inductance of 
the wire per unit length, c is the speed of light, r is 
the wire radius, S is the effective area of the wire cross- 
section, S ~ 7rr^, for (5 > r and S ~ 7r(5(2r — S), for 
S < r, where S = cj ^ (27r crw) is the thickness of the skin 
layer. The average current density in the unit sell can be 
written in the form. 



(2) 



For the waves polarized along the wires the average 
macroscopic electric field E is approximately equal to the 
local field E'. Taking into account the general relation 
between the electric field E and the electric induction D, 
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dielectric permittivity 



(4) 



where ujp ~ {c/d)[2TT/ \n{d/r)]^^'^ is the effective plasma 
frequency, and 7^ — c^/2aShi (d/r). The second term in 
the right-hand side of Eq. Q is in a full agreement with 
the earlier result obtained by Pendry and co-authors P]. 
One should note that the low losses case, i.e. 7^ ^ w, 
corresponds to the condition S <^ r. 

The analysis becomes more involved for calculating a 
nonlinear magnetic response of the composite structure, 
which is determined by the intrinsic properties of the in- 
teracting nonlinear oscillators in the presence of an exter- 
nal periodic force. For the structure under consideration, 
the current induced in each resonator can be found as 



/ = —ir 



- |Ho|Z- 

c - 



(5) 



where Hp is the amplitude of the external field applied to 
SRR, Z = iLuL + R+ {iujC)~^ is the SRR impedance, and 
other parameters are marked in Fig. Q a is the radius, 
L, i?, and C are inductance, resistance and capacitance, 
respectively. The amplitude of the electric field in the 
gap of SRR can be found approximately as follows 



|Eo 



jCda 



(6) 



where dg is the size of the SRR gap. Nonlinear effects 
in Eqs. jSJ, lO appear due to the capacitance C which 
is proportional to e^diSgl^). Therefore, Eq. © gives 
an implicit relation between the amplitude of the local 
electric field in the gap and the amplitude of the (external 
to SRR) magnetic field. 

The relation between the local (microscopic) and av- 
erage (macroscopic) magnetic fields can be obtained by 
averaging a local field over the volume of the unit cell. 



(B) - (H) + ^^H', 



(7) 



where H' is an additional magnetic field induced by 
the alternating external magnetic field Ho in cylindrical 
SRR, which determines magnetization of the composite, 
and F = na^/d'^. 

The average macroscopic field can be also determined 
as follows, 

H = (H) = (1 - F)Ho + F(Ho + H') = Hq + FH'. (8) 

Taking into account that H' = outside SRR, from the 
boundary conditions we obtain the relation 
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(9) 



we can obtain the expression for the effective nonlinear 



where jg is an equivalent surface current in SRR, which 
is equal to the current per unit length. From Eq. |SJl and 



(c) 




0.0001 0.0002 0.0003 0.0004 0.005 0.01 0.015 0.02 

IhP/E ^ IHI^/E ^ 




0.0002 0.0003 
2 2 

IHI/E 



0.01 0.015 

2 2 

IHI TE 



FIG. 2: Real part of the effective magnetic permeability vs. 
intensity of the magnetic field: (a) f2 > f , q = I; (b) f2 < f , 
a = 1, (c) n > 1, a = -1; and (d) f2 < f, a = -I. Black 
- the lossless case (7 = 0), grey-the lossy case (7 = 0.05). 
Dashed curves show unstable branches. 



FIG. 3: Imaginary part of the effective magnetic permeability 
vs. intensity of the magnetic field for 7 = 0.05: (a)f2>l, 
a = 1; (b) n < 1, a = 1, (c) n > 1, a = -1; and (d) Q.<\, 
Of = — 1. Dashed curves show unstable branches. 



Eqs. Q-linj), we obtain an explicit expression for the ef- 
fective magnetic permeability of the composite structure 
(for F < 1): 



/Xeff(H) = 1 



(10) 



where 



W - f)/„d(|e,(H)P)] 

is the eigenfrequency of oscillations in the presence of 
the external field of a finite amplitude, h is the width 
of the ring, F = c^/27r(l — F)aah, for h < S, and 
r = c^/27r(l — F)aa6, for h > S. It is important to 
note that Eq. (|IO|l has a simple physical interpretation: 
The resonant frequency of the artificial magnetic struc- 
ture depends on the amplitude of the external magnetic 
field and, in turn, this leads to the intensity-dependent 
function ficS- 

To be more specific, we consider the Kerr nonlinearity 
of the dielectric in the composite material, i.e 



(11) 



where Ec is a characteristic electric field, a = ±1 stands 
for the focusing or defocussing nonlinearity, respectively. 
Then, the relation between the macroscopic magnetic 
field and the dimensionless nonlinear resonant frequency 
can be obtained from Eqs. ll5l)-(|51), and (|11() as 

|H|2 = a^2£;2 V 1 M ^ (12) 

where = 16e||Q(l — F)uj1h? /c^, O = w/wq, t^o = 
(c/a) [dg/7r(l — F)/ie£)o]^/^ is the eigenfrequency of the 



system of SRRs in the linear limit, X = luonl/^^o, and 
7 — T/uJo- Therefore, we find that the dimensionless 
eigenfrequency of the SRR lattice X is a multi- 

valued function of the magnetic field. This result reveals 
a general property of nonlinear oscillators with a high 
quality factor [Tl| . 

Parametric dependence of the effective magnetic per- 
meability on the magnetic field is determined completely 
by Eqs. H12|) and (|IO|l . Figures |21 and 13 summarize dif- 
ferent types of nonlinear magnetic properties of the com- 
posite, which are defined by the dimensionless frequency 
of the external field il, for both focusing [Fig. IHa,b) and 
Fig. |3Ia,b)] and defocussing [Fig.|2Ic,d) and Fig. |Sl(c,d)] 
nonlinearity of the dielectric. 

In the case oi focusing nonlinearity (i.e. when a — 1), 
the SRR eigenfrequency decreases with the intensity of 
the electromagnetic field because of a growth of the SRR 
capacitance. Then, for fl > 1, the effective magnetic 
permeability of the composite structure grows with the 
field intensity, as shown in Fig.|2Ia). If in the linear limit 
the composite material is left-handed, i.e. i?e(^efi) < 0, 
it will become right-handed for higher intensities of the 
magnetic field [in this reasoning we assume Re{ees) < 0]. 

More complicated behavior of the magnetic permeabil- 
ity is observed for < 1; this is shown in Fig. EJb). 
Here, in the linear limit the real part of the magnetic 
permeability is always positive, but the eigenfrequency 
of SRR decreases with the growth of the magnetic field, 
thus driving the system into the resonance. Since a non- 
linear oscillator has a hysteresis structure of its response 
with a change of an external force (see, e.g., Ref. pdj). 
this leads to multi- valued dependencies. In our prob- 
lem, the nonlinear eigenfrequency is a three-valued func- 
tion of the external magnetic field, and this results in 
jumps of the magnetic permeability with the growth of 
the magnetic field. As follows from Fig. [S^b), the mag- 
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netic field intensity displays a jump of the magnetic per- 
meability from positive to negative values at some -ffci- 
Thus, the initially opaque medium with positive refrac- 
tion becomes a negative-refraction transparent medium 
with the growth of the field intensity. This effect can be 
treated as the second-order phase transition induced by 
the external electromagnetic field. Reversed transition 
takes place when the magnetic field intensity decreases 
to the value Hc2 < Hd- 

In the case of defocussing nonlinearity (i.e. when 
a = —1), the SRR eigenfrequency increases with the am- 
plitude of the external field. That is why the resonance 
effects take place for O > 1, as shown in Figs. I2c,d) 
and Figs. |SIc,d). Here, we observe the opposite behavior 
when the transition from the case i?e(/icff) < to the 
case Re{^cs) > takes place at high values of the exter- 
nal field, and the reversed transition occurs at lower field 
intensities. In the latter case, i?e(/icff) is always positive 
for f7 < 1, see Fig.l^Jd). 

Our results show that the imaginary part of the effec- 
tive magnetic permeability, which determines the struc- 
ture losses, can be controlled rather effectively by a 
proper choice of the intensity of the external high- 
frequency magnetic field, see Fig. |3| We believe that 
this feature may be important for the future applications 
of left-handed materials. 

Due to the high values of the electric field in the gap of 
SRR and the resonant interaction of the electromagnetic 
field with the SRR lattice, the characteristic magnetic 
nonlinearity in such structures is much stronger then the 
corresponding electric nonlinearity. Therefore, the mag- 
netic nonlinearity should dominate in the composite ma- 
terials which display the phenomenon of negative refrac- 
tion. Moreover, nonlinear elements can be used only for 
filling the gap in SRR, allowing an easy tuning by apply- 
ing an external field. 

The possibility of strongly enhanced nonlinearities in 
left-handed metamaterials revealed here may lead to an 



essential revision of the concepts based on the linear the- 
ory, since the electromagnetic waves propagating in such 
materials always have a finite amplitude. At the same 
time, the engineering of highly nonlinear composite ma- 
terials will open a number of novel opportunities for their 
microwave applications, such as frequency multipliers, 
beam spatial spectrum transformers, switchers, limiters, 
etc. 

It is important to note that the hysteresis behavior 
with jumps in the dependencies of the effective material 
parameters has been described above for the stationary 
processes only. Such transitions will display a character- 
istic scale in time or space for initial or boundary prob- 
lems, respectively. Such spatial or temporal scales are 
determined by the relaxation micro-processes in the SRR 
lattice. 

In conclusion, we have presented, for the first time to 
our knowledge, a systematic analysis of nonlinear proper- 
ties of the microstructured materials which display neg- 
ative refraction, the so-called left-handed metamaterials. 
We have shown that the composite metamaterials com- 
posed of a lattice of wires and split-ring resonators pos- 
sess an effective magnetic permeability which depends 
on the intensity of the macroscopic magnetic field in a 
nontrivial way. The magnetic nonlinearity is found to 
be much stronger than the nonlinearity in the dielectric 
properties due to the field enhancement in the split-ring 
resonators. The dependence of the effective magnetic 
permeability on the field intensity allows switching be- 
tween its positive and negative values, i.e. a change of 
the material properties from left- to right-handed and 
back. Such processes can be treated as the second-order 
phase transitions induced by the variation of an external 
electromagnetic field. 
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